Gamma Rays from Kaluza-Klein Dark Matter 
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A TeV gamma-ray signal from the direction of the Galactic center (GC) has been detected by the 
H.E.S.S. experiment. Here, we investigate whether Kaluza-Klein (KK) dark matter annihilations 
near the GC can be the explanation. Including the contributions from internal bremsstrahlung as 
well as subsequent decays of quarks and r leptons, we find a very flat gamma-ray spectrum which 
drops abruptly at the dark matter particle mass. For a KK mass of about 1 TeV, this gives a good 
fit to the H.E.S.S. data below 1 TeV. A similar model, with gauge coupling roughly three times as 
large and a particle mass of about 10 TeV, would give both the correct relic density and a photon 
spectrum that fits the complete range of data. 

PACS numbers: 95.35.+d, 04.50.+h, 98.70.Rz 



INTRODUCTION 

Currently there is great interest in the very high-energy 
gamma-ray signal from the center of the galaxy, where re- 
cently the CANGAROO Q, VERITAS [3 and H.E.S.S. 
3] collaborations have reported results. The three exper- 
iments find different spectra, something that may depend 
on different angular acceptances, energy calibrations, etc. 
The H.E.S.S. result superficially looks more detailed and 
accurate and is therefore the one we will focus on here, 
although of course the three groups will have to find a 
conclusive reason for the reported discrepancies. There 
have already been speculations of a TeV-scale dark mat- 
ter (DM) particle giving, through annihilations, a steady 
source of these photons 0, 0>1H . Of course, there are 
other possibilities for emission of these gamma rays re- 
lated to the physics around the Galactic massive black 
hole 0]. 

The prototype DM particle is the neutralino, the light- 
est supersymmetric particle. The phenomenology of this 
candidate has been worked out in detail over the last two 
decades (for a summary, see Q). However, these mod- 
els give a rather soft gamma-ray spectrum, so one has to 
push masses up to the tens of TeV range to match the re- 
sults of H.E.S.S., which will be our main concern in this 
Letter. We focus instead on an interesting alternative 
which has received much attention in recent years; that 
of universal extra dimensions (UED) IslSLLlJI- I n these 
models, all standard model fields propagate in the higher 
dimensional bulk. For the effective four-dimensional the- 
ory, this means that all particles are accompanied by a 
tower of increasingly more massive Kaluza-Klein (KK) 
states. Momentum conservation in the extra dimensions, 
which corresponds to KK mode number conservation, is 
broken by an orbifold compactification which projects 
out unwanted degrees of freedom at the zero mode level. 
If the boundary terms introduced at the orbifold fixed 
points are identical, then a remnant of KK mode num- 
ber conservation is left in the form of KK parity, and the 
lightest KK particle (LKP) is stable '9]. This is anal- 



ogous to conserved R-parity in supersymmetric models 
which ensures the stability of the lightest supersymmet- 
ric particle. If the LKP is also neutral and non-baryonic, 
it is a potential dark matter candidate. 

We consider the simplest, five-dimensional model with 
one compactificd UED on an S 1 /Z2 orbifold of radius R. 
At tree level, the nth KK mode mass is then given by 
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where ttiew is the zero mode mass. However, the iden- 
tification of the LKP is nontrivial because radiative cor- 
rections to the mass spectrum of the first KK level are 
typically larger than the corresponding electroweak mass 
shifts given by mEw- A one- loop calculation shows 
that the LKP is well approximated by the first KK mode 
of the hypercharge gauge boson B^\ The B^ 1 relic den- 
sity was determined in d0|. Depending on the exact 
form of the mass spectrum and the resulting coannihi- 
lation channels, the limit from the Wilkinson Microwave 
Anisotropy Probe (WMAP) of J] CDM /i 2 = 0.12±0.02 
corresponds to 0.5 TeV < m B (i> < 1 TeV. Here Ocdm is 
the ratio of DM to critical density and h is the Hubble 
constant in units of 100 kms -1 Mpc -1 . Collider mea- 
surements of electroweak observables give a current con- 
straint of R- 1 > 0.3 TeV U EH, whereas LHC should 
probe compactification radii up to 1.5 TeV [l3| . 

The prospects of KK DM detection have been studied 
in some detail 0, 0, 0, 0, 0, . Secondary gamma 
rays from cascading quark decays due to B^'B^- 1 ' annihi- 
lation in the GC produce relatively soft spectra fldl llij i . 
However, unlike the supersymmetric case, charged lepton 
production is not helicity suppressed, which, e.g., results 
in a striking peak signal in the positron spectrum 0,0] 
for masses < 0.5 TeV. This has led us to also study the 
decays of r leptons and the spectrum of primary contin- 
uum gamma rays from the production and subsequent 
radiative photon emission of high-energy charged leptons. 
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PRIMARY AND SECONDARY GAMMA RAYS 
FROM B W B W ANNIHILATIONS 



We first consider primary gamma rays. At tree level, 
with all other first level KK modes degenerate in mass, 
pairs annihilate into quark pairs (35%), charged 
lepton pairs (59%), neutrinos (4%), charged (1%) and 
neutral (0.5%) gauge bosons, and Higgs bosons (0.5%). 
These branching ratios agree with those obtained by 
Servant and Tait |l0l | (however, they do not take into 
account electroweak symmetry breaking). In this Let- 
ter we only consider photons radiated from charged lep- 
tons We do not concern ourselves with the radia- 
tive processes which electrons and positrons encounter 
during their propagation through the Galaxy, as these 
give mostly low-energy photons. The exception may 
be inverse Compton scattering on infrared photons and 
starlight, b ut g enerally this is expected to be only a small 
correction [20j . 




The tree level Feynman diagrams which contribute to 
the process B^'B^ — » £ + £~j are shown in Fig.^ The 
computation of the cross section is straightforward and 
we find that the differential photon multiplicity is well 
approximated by 
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where x = E 7 /m B d) . From the electromagnetic coupling 
and the phase space difference between two- and three- 
body final states, one would expect an answer of the or- 
der of a/ir times a large logarithm, which is related to 
a collinear divergence. Our calculation shows that there 
is indeed such a leading logarithmic term, giving large 
contributions for high photon energies £L. Restricting 
ourselves to these energies (at lower energy, quark frag- 
mentation dominates anyway), the radiative correction, 
although large, is not more than some five percent of the 
lowest order cross section. Therefore, there is no need in 
this first treatment of the problem to sum leading loga- 
rithms, but this could eventually be done [2l| . 

We would like to emphasize that our result (0 is al- 
most entirely due to the very large mass of the B^ 
and practically independent of the initial state spin. Ex- 
pressed in the scaling variable x, it is furthermore quite 
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FIG. 2: The total number of photons per B*- 1 ' annihila- 
tion (solid line), multiplied by x 2 = (2J 7 /m s (i) ) 2 . Also shown 
is what quark fragmentation alone would give (dashed line), 
and adding to that the contribution from r leptons (dotted 
line). Here we have assumed a B^ mass of 0.8 TeV and a 
5% mass splitting at the first KK level, but the result is quite 
insensitive to these parameters. 



insensitive to the B^ 1 ' mass and the mass splitting at the 
first KK level. This result therefore applies not only to 
KK DM, but to any weakly interacting massive parti- 
cle (WIMP) candidate with nonsuppressed couplings to 
charged leptons. 

Let us now consider secondary gamma rays. As already 
mentioned, cascading decays of qq final states have been 
studied previously [lj, |l8j. Here, we also include the 
semihadronic decays of r leptons, which are important in 
KK DM models since they have a fairly hard spectrum 
and a branching ratio of around 20%. We will use the re- 
cent results of Fornengo, Pieri, and Scopel [22|, who have 
used the Pythia Monte Carlo code |23j to parametrize 
dN^ T /dx for quarks and r leptons with a center of mass 
energy of 1 TeV. We neglect the few percent going into 
W, Z, and Higgs final states. 

The total number of photons per B^B^ annihilation 
is given by dN° s /dx = KidN^/dx, where the sum is 
over all processes that contribute to primary or secondary 
gamma rays, and k% are the corresponding branching ra- 
tios. The result is shown in Fig. [21 Previous analyses 
of the photon flux correspond to the much softer and 
sharply falling spectrum from quark fragmentation alone. 
Also shown is the more important contribution from r de- 
cays, and of course the radiative direct process which is 
the main topic of this Letter. 



GAMMA-RAY FLUX FROM THE 
GALACTIC CENTER 



The details of the Galactic halo profile are to a large 
extent unknown. High resolution N-body simulations fa- 
vor cuspy halos, with radial density distributions ranging 
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FIG. 3: The H.E.S.S. data 0] compared to the gamma-ray 
flux from a region of 10 -5 sr encompassing the GC, for a B^ 1 ' 
mass of 0.8 TeV, a 5% mass splitting at the first KK level, 
and a boost factor b around 200 (dashed line) . The solid line 
corresponds to a hypothetical 10 TeV WIMP with similar 
couplings, a total annihilation rate given by the WMAP relic 
density bound, and a boost factor around 1000. 



from r _1 [24j to r -15 [25| in the inner regions. It has 
further been suggested that adiabatic accretion of DM 
onto the massive black hole in the center of the Milky 
Way may produce a dense spike of r~ 2A 26]. This has, 
however, been contested [23, |28|. On the other hand, adi- 
abatic contraction from the dense stellar cluster, which is 
measured to exist near the GC, is likely to compress the 
halo DM substantially [2^, I^O] • This means that there 
is support for a rather dense halo profile very near the 
center - something that may be tested with the new 30- 
m-class telescopes [Ul. Bearing these uncertainties in 
mind, we will use the moderately cuspy profile by 

Navarro, Frenk and White (NFW) jH]. 

Following [32l |. the differential gamma-ray flux from 
WIMP annihilation in the GC can be written as 
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where <7tot^ is the total B^ annihilation rate, and 
(^Gc)Afi is a dimensionless line-of-sight integral aver- 
aged over Af2, the angular acceptance of the detector. 
For a NFW profile, one expects (JGc)AnAtt = 0.13 b for 
Ail, = 10" 5 0]» which is appropriate for the H.E.S.S. 
telescope. Here, we allow for an explicit boost factor b 
that parametrizes deviations from a pure NFW profile 
(b = 1) and may be as high as 1000 when taking into ac- 
count the expected effects of adiabatic compression [3(J ■ 
For KK DM, a tot v ~ 3 • 10- 26 (0.8 TeV/m B( i)) 2 cm 3 s" 1 

M 

In Fig. |3 we plot the total differential gamma-ray 
flux for a mass of 0.8 TeV together with the re- 



cent H.E.S.S. observations with its surprising, relatively 
hard spectrum As can be seen, an excellent fit to 
the data is obtained for energies up to m B (i), requiring 
only a moderate value for the boost factor b of roughly 
100. A more detailed analysis would of course take the 
energy resolution of the detector into account, giving a 
less rapid cutoff at the B^ mass. 

In order to explain the whole observed spectrum with 
KK DM, one would thus need a mass of m B (i) > 7 TeV. 
Such a high mass, however, is outside the range selected 
by the relic density constraint and thus disfavored. It 
has to be kept in mind, though, that models with ex- 
tra dimensions may have quite a different behavior near 
freeze-out than ordinary models, which could alter the 
conclusions about the preferred mass range. Abnormally 
large running of the U(l) coupling might also allow for 
larger masses 0- 

As remarked before, we expect a similar spectrum for 
any WIMP dark matter candidate with nonsuppressed 
couplings to charged leptons. Such a model would give 
a nice fit to the full H.E.S.S. data, if the canonical 
WIMP relation given by the WMAP relic density bound, 
(ctw)wimp ~ 3 • 10~ 26 cm 3 s _1 0, allows for particle 
masses above 7 TeV. For the KK DM scenario consid- 
ered here, this would require the B^'s gauge coupling 
to be increased by a factor of roughly three. As an illus- 
tration, we have included such a hypothetical particle of 
mass 10 TeV in Fig.|3with a boost factor of roughly 1000. 
Models which only get a contribution from quark frag- 
mentation, on the other hand, need a much larger boost 
factor to fit the observed range. For supersymmetry 
for example, one is furthermore pushed to a neutralino 
mass larger than about 12 TeV, which starts to involve 
severe fine-tuning of parameters. It has to be remem- 
bered, though, that the radiative effects investigated in 
this Letter have not yet been systematically studied for 
supersymmetric models |2l| . 

Of course, the models presented in Fig. |31 must not vi- 
olate the constraints given by the existing EGRET j3f| 
measurements in the energy range 30 MeV - 10 GeV. 
We find, however, that these upper bounds are easily 
satisfied. There are some further tests which should be 
performed in the future to check the KK hypothesis of 
the TeV gamma-ray signal. Of course, a signal from an- 
nihilating DM is not expected to show any detectable 
time dependence. Since the mass of the DM particle is 
an abolute cutoff for the energy of any secondary parti- 
cle, a steplike drop should definitely be searched for. If 
one gets instruments with better energy resolution, the 
line from B^B^ —> 77 annihilation \3(i\ may form the 
final proof of the particle DM hypothesis . It will be 
interesting to follow the wealth of new results expected 
from H.E.S.S., as well as MAGIC 0, CANGAROO Q 
and VERITAS 0]. The GLAST satellite [H should also 
give important data in the range where the various fi- 
nal states cooperate to give a signal. The fate of the 
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produced electrons and positrons should also be investi- 
gated, something that will need educated guesses for the 
magnetic field and exact DM density distribution near 
the GC. In fact, a slightly boosted NFW profile seems 
to predict a radio signal which is quite close to that ob- 
served 0, |2(J • We find further that the direct anni- 
hilation into neutrinos will, in the KK model and with 
the same halo parameters explaining the gamma-ray sig- 
nal, give an event rate corresponding to a few events per 
km 2 per year ^^|, something that could eventually be 
detected in a neutrino telescope with good view of the 
GC. 



CONCLUSIONS 

In this Letter we have presented the somewhat surpris- 
ing result that gamma rays radiated from charged leptons 
in B^B 1 - 1 ' annihilations (e.g., near the GC, but also in 
other overdense regions of the halo, and maybe even in 
the diffuse extragalactic gamma-ray spectrum 39]) could 
be the most promising way to detect the effects of KK 
DM and to possibly differentiate it from other candidates. 
We have seen that the signature is a surprisingly hard 
spectrum, at high energies several orders of magnitude 
larger than the previously studied case of quark frag- 
mentation alone, which extends all the way to the mass 
of the DM particle and then drops rapidly. Even if the 
present gamma-ray signal detected by H.E.S.S. will turn 
out to have another cause, the sharp drop at E 7 = m B (i) 
is a nice signature for that and other experiments to look 
for. 

L.B. is grateful to the Swedish Science Research Coun- 
cil (VR) for support. We acknowledge useful discussions 
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